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Purpose: To improve adaptive reconstruction of multichannel MR images by simultaneously remov- 
ing nonsmooth phase and signal-loss imaging artifacts. 

Methods: The improved adaptive reconstruction consists of three steps: (1) modified multichannel 
images are first derived by dividing raw multichannel images by a reference image (i.e., a normal- 
ized single-channel image); (2) the modified multichannel images are smoothed by a low-pass filter; 
(3) adaptive spatial matched filters determined from the smoothed multichannel images are utilized 
to obtain multichannel combined images. Numerical simulations, as well as MRI experiments, on 
phantoms and human subjects are performed to evaluate and compare the effectiveness of this improved 
adaptive reconstruction approach against traditional coil combination methods. 

Results: Both simulation and MRI experimental results demonstrated that the proposed improved 
adaptive reconstruction method is able to obtain combined images with reduced nonsmooth phase 
and signal-loss imaging artifacts. 

Conclusions: A novel multichannel image reconstruction method is developed that produces high 
quality multichannel combined images. © 2015 American Association of Physicists in Medicine. 


[http://dx.doi.org/10.1118/1.4905163] 
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1. INTRODUCTION 


Multichannel receiver coils have been widely used for MRI 
since 1990s.! Due to the inhomogeneous coil sensitivity map 
of surface coils, the object signal intensity rapidly decreases 
with increased distance from the coil location. Thus, a proper 
multichannel image reconstruction method is needed to 
generate a combined image with a uniform intensity. 

For multichannel magnitude image reconstruction, the 
commonly practiced method is the sum-of-squares (SOS) 
approach.' However, the SOS method is known to reduce 
imaging signal-to-noise ratio (SNR) and introduce signal bias: 
signal appears larger than it should be, especially in regions 
with low SNR.” To address these problems, various coil sensi- 
tivity based methods for multichannel image combination 
have been developed, such as electromagnetic field calcula- 
tion methods,’ reference scan methods,*> and multichannel 
images based methods.°” On the other hand, for multichannel 
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phase image reconstruction, both the weighted mean (WM) 
methods*!! and the coil sensitivity based methods**’ are 
commonly used. The SNR of the phase images produced 
by the coil sensitivity based methods are higher than that 
of the WM methods, especially when the SNR is low.’ This 
work focuses on improving one of the coil sensitivity based 
methods, the adaptive reconstruction (AR) method proposed 
by Walsh et al.,’ for both the multichannel magnitude and 
phase image reconstruction. AR is the default multichannel 
image combination method used in Siemens scanners. 

The AR method uses a statistical description of the ex- 
pected signal and noise covariances to derive estimates of 
the optimal spatial matched filters (SMF, i.e., coil sensitivity 
weights). The SNR of both combined magnitude and phase 
images produced by AR method is improved compared to 
conventional SOS and WM methods, respectively. However, 
SMF of AR tend to result in inconsistencies and sometimes 
failure in phase determination (nonsmooth phase artifacts), 
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particularly at points where the phases between the channels 
vary substantially to cause signal cancelation. !? In practice, 
nonsmooth phase artifacts often occur in the combined image 
when multichannel images are in regions with low SNR. 

Walsh et al.’ mentioned that SMF computed from low-pass 
filtered multichannel images may eliminate/reduce their noise 
correlation with raw multichannel images. This modified AR 
(MAR) method could further reduce background noise in the 
combined magnitude image and increase the SNR of the com- 
bined phase image. Besides, the nonsmooth phase artifacts in 
the combined image could be significantly reduced because of 
the increased phase consistencies of the low resolution SMF 
with improved SNR. However, in practice, low-pass filtered 
multichannel images may lose regional coil sensitivity infor- 
mation when the underlying regional MR signals suffer rapid 
phase variation due to the strong susceptibility inhomogeneity, 
such as in the air—tissue interfaces.’-!3!4 Then, the combined 
image using SMF determined from the filtered multichannel 
images will suffer signal-loss artifacts both in magnitude and 
phase images. 

In terms of physical nature, the phase map of each multi- 
channel image represents mixtures of RF offset phase and 
inhomogeneous susceptibility components.'? Except the RF 
offset phases induced by the different receive coils, the other 
phase components including static field inhomogeneity (due 
to imperfect shimming) and time varying inhomogeneity (due 
to eddy currents and object motion) are exactly the same in 
all the multichannel images. In this note, we develop a robust 
modified adaptive reconstruction (ROMAR) method to elim- 
inate the signal-loss artifacts in the MAR combined image 
due to the existence of large inhomogeneous susceptibility 
components in the raw multichannel phase images. The key 
steps of ROMAR are as follows: (1) the large inhomogeneous 
susceptibility components in the raw multichannel phase im- 
ages are removed by dividing raw multichannel images by a 
reference image (i.e., a normalized single-channel image); (2) 
the obtained new multichannel images without rapid phase 
variation patterns are smoothed by a low-pass filter, and the 
filtered multichannel images contain sufficient coil sensitivity 
information for adaptive constructing the SMF, which have 
a low noise correlation with raw multichannel images; (3) 
both the magnitude and phase components of the combined 
image using the above SMF are SNR-optimal and free of both 
nonsmooth phase and signal-loss imaging artifacts. 


2. THEORY 


Complex channel image C;(r) acquired by j-th receive coil 
can be described as a pixel by pixel product of native MR 
signal p(r) multiplied by corresponding coil sensitivity S;(r) 
with additive noise n,(r), 


Cir) = pr) Sir) +n,(r), 


where N is the total number of receive channels and r is the 
pixel location. The combined image is the linear superposition 
of weighted multichannel images! 


I(r) =W(r)"C(r), (2) 


j=1,2,..., N, (1) 
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where /(r) (r is indexing the pixel) is the combined image. 
W(r) is a vector of N channel components and vector cell is the 
coil sensitivity weight W,(r). H denotes the Hermitian trans- 
pose operation. C(r) is also a vector of N channel components 
and vector cell is the channel image C;(r). 

AR method demonstrates that the optimal spatial matched 
filter is obtained by maximizing imaging SNR, and the SNR 
formulation is expressed as follows:’ 


SNR = ——— (3) 


where R, and R,, are the second order array correlation matrices 
of the MR signal and noise, respectively. The matrix cells of Rs 
and R, are expressed as follows:’ 


RZK)= E (CAr)Cry), fal, N, k=1,..., N, (4) 
(reSROD 

RGk)= E (nry), fal... N, k=1,..., N, (5) 
(reSROI) 


where SROI is a specified set of pixel coordinates and * repre- 
sents the Hermitian conjugate operator. Operator E(-) is the 
signal’s expectation. Then, the optimal spatial matched filter 
W(r) can be generated from the eigenvector corresponding to 
the maximum eigenvalue of the matrix product 


Pak R;: (6) 


The noise correlation matrices R„ can be derived from a pres- 
can noise calibration or can be simply assumed to be an 
identity matrix. Then, the obtained W(r) is used pixel by pixel 
or block by block basis to combine multichannel images.’ 

Further SNR improvement can be obtained by MAR which 
utilizes smoothed (i.e., low-pass filtered) multichannel images 
to compute SMF.’ Commonly, low-pass filtering is applied in 
the image domain by windowing the k-space data.'>!® The 
native MR signal (p(r)) may have both low and high spatial 
frequency components (p(r)!"4 and p(r)hiz'fea) and coil 
sensitivity profile S(r) only has a low frequency character. 
Then, both original and smoothed multichannel images can be 
expressed as follows: 


Cir) = plr) ESC) + plr) eaS) +0,(r), 
j=1,.., N, (7) 
(CE) = pr 48 (r) + (nr), fH... N, (8) 
where ( ) indicates that the image has been smoothed by a low- 
pass filter and it is not the signal’s expectation. If the native 
MR signal has enough low frequency components (p(r)'°¥ #4) 
in each pixel location, coil sensitivity information is still re- 


mained in (C,(r)). The corresponding signal and noise array 
correlation matrix R, and R, can be expressed as 


Rs(j,k)= E  (Cj(r)XCk(r))), 
(r €SROI) 


j=1,.., N,k=1,..., N, (9) 
RG, = E (aY), 
(r e(SROD) 
j=1,.., N, k=1,..., N. (10) 


Most noise energy in original multichannel images is re- 
moved by the smoothing operation. Thus, the derived SMF 
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Fic. 1. Reconstruction results of noisy multichannel phantom images. The SOS image is shown in the first row. The second to fourth rows show the magnitude 
and phase results of AR, MAR, and ROMAR methods, respectively. Images in the first to fourth columns are the magnitude, difference, phase, and high-pass 
filtered phase maps of different combination methods, respectively. High-pass filtered phase maps are shown in the range of [—1,1] (the unit is radian). Imaging 


artifacts are shown in the red circles in the above images. 


determined from matrix R, exhibit almost no correlation with 
the original multichannel images in the nonsignal regions. So 
the combined image has a further reduced noise floor. 


2.A. ROMAR of multichannel MR images 


The image combined by MAR method always produces 
an optimized SNR. However, some regional MR signals may 
exhibit rapid phase variation due to the strong susceptibility 
inhomogeneity, and these regions almost do not contain the 
low frequency native MR signal components. Thus, the low- 
pass filtered coil image will lose coil sensitivity information 
in these rapid phase variation regions and consequently the 
combined image will suffer signal-loss artifacts. 

In light of the problem in MAR, ROMAR is introduced 
to avoid the signal-loss artifacts and simultaneously optimize 
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SNR of the combined image. ROMAR removes the large inho- 
mogeneous susceptibility components in raw multichannel 
phase images through dividing multichannel images by a 
reference image. Then, the new channel image C Rr) does not 
contain rapid phase variation patterns and it is expressed as 
follows: 


CR(r) =Ci(r)/normalize(Cr(r)), i=1,...,N, RELI, N], 
(11) 


where the reference image (normalize(Cr(r))) can be derived 
by a normalized single-channel image Cr(r) (i.e., normalize 
(Cr(r)) = exp(ig(Cr(r))), where y(Cr(r)) is the phase of 
Cr(r)). The phase components remained in the new chan- 
nel image CR(r) are only RF offset phase induced by the 
receive coil, and the RF offset phases are always globally 
smooth. Then, C? (r) is smoothed by a low-pass filter. The coil 
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Fic. 2. Reconstruction results of noisy multichannel head images. The SOS image is shown in the first row. The second to fourth rows show the magnitude and 
phase results of AR, MAR, and ROMAR methods, respectively. Images in the first to fourth columns are the magnitude, difference, phase, and high-pass filtered 
phase maps of different combination methods, respectively. High-pass filtered phase maps are shown in the range of [—1,1] (the unit is radian). Imaging artifacts 


are shown in the red circles in the above images. 


sensitivity information is well preserved in the filtered image 
(CX(r)). The correlation matrix Rs is similarly calculated by 
Eq. (9). In addition, a constant divisor (normalize(Cr(r))) 
does not change the physical content of the matrix cells 
of R;. Finally, the SMF determined from R, are used for 
multichannel image combination. For ROMAR, since the 
smoothing procedure has been employed to the multichannel 
images without rapid phase variation, signal-loss artifacts are 
effectively avoided. 


3. MATERIALS AND METHODS 


To validate the proposed method, both phantom and in vivo 
experiments were performed. The multichannel images were 
combined, and the results were compared to those of SOS, AR, 
and MAR techniques. In addition, the proposed method is also 
compared with the WM method in phase reconstruction. 
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3.A. MR experiments 


A group of ACR phantom and human brain data were 
acquired on a Siemens 3T scanner using a 12-channel head 
coil with a 2D gradient echo sequence (TR/TE/@ = 620 ms/ 
20 ms/3°, FOV = 220 x 220 mm?, matrix = 256 x 256, slice 
thickness = 3 mm, slice number = 18) to produce noisy data. 
The low SNR data were used to test AR, MAR, and ROMAR 


Taste I. The noise level of the top left regions (size: 40X40) in both phantom 
and head images. 


sos AR MAR ROMAR 

Piatt Mean values 0.153 0.088 0.075 0.063 
antom SD values 0.041 0.046 0.040 0.032 

Hed Mean values 0.165 0.096 0.081 0.076 
a SD values 0.044 0.050 0.042 0.036 
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Fic. 3. Reconstruction results of multichannel 3D SPGR images by MAR (first and third rows) and ROMAR (second and fourth rows) methods. Images in the 
first to fourth columns show the magnitude, difference, phase, and high-pass filtered phase maps, respectively. The images shown in the first and second two 
rows are two different slices acquired with same sequence. High-pass filtered phase maps are shown in the range of [—1,1] (the unit is radian). Imaging artifacts 


are shown in the red circles in the above images. 


methods. The second head data were acquired with a 3D full 
flow-compensated SPGR sequence (TR/TE/@ = 40 ms/32 
ms/15°, FOV = 230x230 40 mm, matrix = 256 x 256 x 40). 
The high SNR SPGR data were used to compare MAR and 
ROMAR methods. Written informed consent was obtained 
from all subjects prior to the examination. In addition, human 
brain data in low SNR were also acquired with a 2D dual echo 
SPGR sequence (TR/TE,/TE2/@ = 40 ms/21.8 ms/24.9 ms/5°, 
FOV = 240 x 240 mm?, matrix = 256 x 256, slice thickness 
=2 mm, slice number = 18) to compare the proposed method 
with the WM method in phase difference reconstruction.*? 


3.B. Data processing 


All methods were implemented in MatLaB 2013a (64-bit; 
Mathworks, Natick, MA). In the data processing of both phan- 
tom and in vivo data, SMF of AR was applied by a pixelwise 
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combination with a size of 16 x 16 SROI correlation matrices. 
The smoothing procedure of MAR and ROMAR was carried 
out by k-space truncation followed by applying a Hanning 
window and zero-filling. Truncation factor (TF) was defined 
by dividing original encoding steps by the remained encoding 
steps after truncation. All the data dimensions had the same 
TF. TFs of both MAR and ROMAR methods were typically 
8.1517 An 8x8 SROI correlation matrix was employed for 
SME calculation of both MAR and ROMAR. 

For ROMAR, any given single-channel image (choose 
channel #1 here) of both phantom and in vivo data processed 
in this note can be employed as the reference image. Though 
SOS is not a SNR optimized coil combine method, it can still 
serve as a useful guide for object signal detection. Difference 
maps were calculated to investigate signal difference between 
SOS images and magnitude components of AR, MAR, and 
ROMAR combined images. High-pass filtered phase images 
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(using typical Gaussian window with o = 42) of AR, MAR, 
and ROMAR results were also computed to identify the 
phase image quality.!* For 3D SPGR data, susceptibility maps 
are determined from the above combined phase images by 
recently developed quantitative susceptibility mapping (QSM) 
method to test the combined phase quality. !° 

For phase difference reconstruction, two groups of the 
multichannel images (the dual echo data) with different echo 
times were reconstructed, respectively, by the proposed RO- 
MAR method. Then, the phase difference image was obtained 
by Eq. (4) of Ref. 9. The WM method using Eq. (3) of Ref. 9 
was employed for comparison. 


4. RESULTS 


The reconstruction results of noisy phantom and head data 
were shown in Figs. | and 2, respectively. The mean and stan- 
dard deviations (SDs) values of the top left noise regions (size: 
40 x40) in magnitude maps of SOS, AR, MAR, and ROMAR 
were shown in Table I. Both the mean and SD values of the 
noise regions in ROMAR magnitude maps were minimum 
when comparing with those obtained from the other methods. 
In addition, the SDs of AR signal were larger than those results 
of SOS, MAR, and ROMAR methods. This occurs because of 
the unstable correlation matrix averaging procedure of AR. As 
can be seen in Figs. | and 2, both phase and high-pass phase 


MAR 


images by AR suffered some nonsmooth phase artifacts (seen 
in the circle regions). The results of MAR incurred in some 
signal-loss artifacts, both in magnitude and phase images, due 
to the rapid phase variation of the data. In contrast, both the 
magnitude and phase images of ROMAR combined image 
were free of artifacts and showed high image quality. 

MAR combined images of 3D SPGR data still presented 
signal-loss artifacts (see the circle regions in Fig. 3). Figure 4 
shows the susceptibility maps which correspond to the com- 
bined phase images in Fig. 3. The small artifacts in the MAR 
phase images caused extremely large errors in susceptibility 
maps determination. 

As can be seen from the results of both noisy phantom and 
head data, nonsmooth and signal-loss artifacts were appeared 
in AR and MAR combined images, respectively. Besides, 
signal-loss artifacts still damaged the MAR combined image 
of high SNR SPGR data. In contrast, whatever images re- 
constructed from noisy or high SNR data, ROMAR method 
always showed SNR-optimized and artifacts-freed combined 
magnitude and phase images. 

The results of phase difference reconstruction were shown 
in Fig. 5. It is seen that signal variance within the circular 
region of the WM phase image is larger than that of the 
ROMAR phase image, which shows that the phase difference 
image reconstructed by ROMAR possesses a higher SNR than 
that of the WM method. 


ROMAR 


Fic. 4. Reconstruction results of susceptibility maps. The two susceptibility maps in the first column are determined from corresponding MAR combined phase 
images in Fig. 3. The susceptibility results of ROMAR are shown in the second column. 
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WM 


ROMAR 


Fic. 5. Reconstruction results of phase difference imaging. The left and right images are phase difference images reconstructed by the WM and proposed 
ROMAR methods, respectively. The signal variance in the red circle regions of the WM and ROMAR phase images are 0.57 and 0.45, respectively. 


5. DISCUSSION 


Both phantom and in vivo studies have demonstrated that 
our proposed image combination method, relative to conven- 
tional image reconstruction methods, yields more accurate and 
robust performance, in particular when multichannel images 
are noisy or oscillating. The developed ROMAR method al- 
ways showed SNR-optimized and artifacts-freed combined 
magnitude and phase images. 

SME of the proposed ROMAR were applied by pixelwise 
combination with 8x8 SROI correlation matrices. Signal aver- 
aging of pixels in SROI box is beneficial for the continuity of 
combined phase image. Similar to the original AR method, the 
size of SROI can be selected adaptively. However, a small size 
of SROI should be sufficient because the majority of the noise 
has been removed by the smoothing method. In our tests, TF 
values from 4 to 20 can be used in practice. A large TF value 
(>28) is recommended for full elimination of noise energy 
and increased phase consistency between channel images. 
Besides, because of the elimination of susceptibility compo- 
nents from raw multichannel phase images, the new channel 
images contain sufficient low frequency components. Thus, 
coil sensitivity information is well preserved in the filtered 
multichannel images [seen in Eq. (8)] even with a large TF. 

The data processed by the proposed ROMAR method in 
this note are fully sampled, which have a higher SNR, and 
do not contain aliasing artifacts compared with undersam- 
pled data. In practice, when the acquired data have a low 
SNR such as diffusion weighted imaging (DWI) with a high 
b-value, or when the imaging contrast is sensitive to both 
aliasing artifacts and noise level such as functional MRI, or 
when the high resolution images which contain small tissue 
details are needed such as whole brain 7; weighted imaging 
for voxel-based morphometric (VBM) analyzing or QSM, we 
recommend these data be fully sampled and combined by the 
proposed ROMAR method to obtain images with high quality. 

The ROMAR method can also be combined with parallel 
imaging techniques such as SENSE (Refs. 4, 15, and 20) and 
GRAPPA.2!:? For combination with SENSE, the SMF (1.e., 
coil sensitivity weights) of ROMAR is derived from the fully 
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sampled center k-space lines.” Then, Eq. (16) of Ref. 20 may 
be employed for optimal magnitude and phase reconstruc- 
tion. To combine with GRAPPA, undersampled data should 
be reconstructed by GRAPPA algorithm first to derive the 
full resolution multichannel data. Then, optimal multichannel 
image combination can be derived by the proposed ROMAR 
method. 

In addition, the proposed ROMAR method inherits the 
capability of suppressing motion-related artifacts in the orig- 
inal AR implementation. Besides QSM, the developed RO- 
MAR method is potentially useful for phase contrast (PC) flow 
imaging. The image reconstruction procedure for PC flow 
imaging is the same as the phase difference reconstruction 
described above. Hence, it is advantageous to use the proposed 
ROMAR method than the WM method? in PC flow imaging, 
especially when the SNR is low. 


6. CONCLUSION 


An improved adaptive reconstruction method is proposed 
in this work to obtain accurate coil sensitivity estimations 
for SNR-optimized and artifacts-freed multichannel image 
combination. Compared to the conventional AR, MAR, and 
WM methods, more robust and reliable results can be obtained 
from the proposed ROMAR method. 
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